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Abstract 

Hyper uniform many-particle distributions possess a local number variance that grows more 
slowly than the volume of an observation window, implying that the local density is effectively ho- 
mogeneous beyond a few characteristic length scales. Previous work on maximally random strictly 
jammed sphere packings in three dimensions has shown that these systems are hyperuniform and 
possess unusual quasi-long-range pair correlations decaying as r -4 , resulting in anomalous logarith- 
mic growth in the number variance. However, recent work on maximally random jammed sphere 
packings with a size distribution has suggested that such quasi-long-range correlations and hyper- 
uniformity are not universal among jammed hard-particle systems. In this paper we show that such 
systems are indeed hyperuniform with signature quasi-long-range correlations by characterizing the 
more general local- volume- fraction fluctuations. We argue that the regularity of the void space in- 
duced by the constraints of saturation and strict jamming overcomes the local inhomogeneity of 
the disk centers to induce hyper uniformity in the medium with a linear small- wavenumber non- 
analytic behavior in the spectral density, resulting in quasi-long-range spatial correlations scaling 
with r~( d+1 ^ in d Euclidean space dimensions. A numerical and analytical analysis of the pore-size 
distribution for a binary MRJ system in addition to a local characterization of the n-particle loops 
governing the void space surrounding the inclusions is presented in support of our argument. This 
paper is the first part of a series of two papers considering the relationships among hyperuniformity, 
jamming, and regularity of the void space in hard-particle packings. 
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I. INTRODUCTION 



Maximally random jammed (MRJ) packings of hard particles are the most disordered 
structures, according to some well-defined order metrics, that are rigorously incompressible 
and nonshearable l|. These systems are prototypical "glassy" structures in the sense that 
they are structurally rigid yet lack Bragg peaks in their scattering spectra j^j. In this sense, 
the idea of the MRJ state has replaced the mathematically ill-defined notion of random close 
packing jl]. Nearly half-a-century ago these systems were thought to describe the disordered 
structure of liquids jsj], but it is now known that three-dimensional MRJ monodisperse 
sphere packings possess unusual quasi-long-range (QLR) pair correlations decaying as r -4 
J]. This property is markedly different from typical liquids, in which pair correlations decay 



exponentially fast {2, 5, 6]. Similar QLR behavior has also been observed in noninteracting 
spin-polarized fermionic ground states [2,8], the ground state of liquid helium 9], and the 
Harrison-Zeldovich power spectrum of the density fluctuations of the early Universe 10]. 
However, for each of these examples and for MRJ hard sphere packings, the structural 
origins of these correlations have been heretofore unknown, even for monodisperse systems. 
Furthermore, it is an open problem to generalize these QLR correlations for MRJ states 
to polydisperse packings, in which the jamming properties are intimately related to the 
size-distribution of the particles [11 ]. 

Motivated by the observation that MRJ packings are structurally rigid with a well-defined 
contact network, Torquato and Stillinger conjectured jij] that all strictly jammed (i.e., me- 
chanically rigid), saturated 12J packings of monodisperse spheres in d- dimensional Euclidean 
space M. d are hyperuniform, meaning that infinite-wavelength local density fluctuations van- 
ish {5}], a proposition for which no counterexample has been found to date 13]. This conjec- 
ture suggests that saturation and strict jamming are sufficient to induce hyper uniformity, 
albeit not necessary 



re 



5|. Hyperuniform systems play an integral role in understanding the 
ationship between fluctuations in local material properties and microstructural order [4- 
16j. These systems have applications to the large-scale structure of the Universe [10] . 
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the structure and collective motion of grains in vibrated granular media [17|, the structure 
of living cells 18], transport through composites and porous media [lil, the study of noise 



and granularity of photographic images [20|, |2l|, identifying properties of organic coatings 
22 1, and the fracture of composite materials [23]. For microstructures consisting of "point" 



particles, one considers fluctuations in the local number density within some observation 
window. Hyperuniform point patterns possess local density fluctuations that asymptotically 
grow more slowly than the volume of the window. Recently, the concept of hyperuniformity 
has been extended to include systems composed of finite- volume inclusions of arbitrary ge- 



ometries in these cases, the quantity of interest is the fluctuation in the so-called local 
volume fraction, defined as the fraction of the volume within an observation window covered 
by a given phase. Hyperuniform heterogeneous media possess local-volume-fraction fluctua- 
tions that asymptotically decay faster than the volume of the observation window, implying 
that the local volume fraction approaches a global value beyond relatively few characteristic 
length scales. 

Previous work on three-dimensional (3D) MRJ monodispere sphere packings j|| has sup- 
ported the Torquato-Stillinger conjecture by showing that the structure factor S(k), pro- 
portional to the scattering intensity, approaches zero linearly as the wavenumber k — > 0, 
inducing a QLR power- law tail r -4 in the pair correlation function g^O")- This behavior im- 
plies that local-number-density fluctuations grow logarithmically faster than the surface area 
of an observation window but still slower than the window volume. In this sense, the "de- 
gree" of hyperuniformity in MRJ pac king s is minimal among all strictly jammed saturated 



packings. However, recent numerical 24J and experimental 25| work on polydisperse MRJ 
packings has suggested that hyperuniform quasi-long-range correlations are not a universal 
signature of the MRJ state. Unlike monodisperse systems, the distribution of particle sizes in 
a polydisperse packing introduces locally inhomogeneous regions as the particles distribute 
themselves through space as shown in Figure [TJ One can see in the binary packing that local 
clusters of small particles are distributed near and around larger inclusions, and the result 
is that the point pattern generated by the disk centers possesses local inhomogeneities that 
are expected to (and indeed do) induce volume-order scaling within the number variance. 
The situation is apparently even more complex for the polydisperse system in Figure [1] since 
the size distribution results in a highly inhomogeneous local structure with small particles 
trapped between larger ones with a high probability. 

These observations have raised a number of quantitatively and conceptually difficult 
questions. First, what is the appropriate extension of the Torquato-Stillinger conjecture 
for monodisperse MRJ packings to systems with a size distribution? Clearly, one must 
explicitly account for the shape information of the particles. Second, in the event that one 



(a) (b) 

FIG. 1. (Color online) (a) A binary packing of hard disks near the MRJ state, (b) A polydisperse 
packing of hard disks near the MRJ state. 



can generalize the Torquato-Stillinger conjecture, there is to date no satisfactory structural 
explanation for the linear small-wavenumber scaling of the structure factor observed for 
3D MRJ monodisperse hard sphere packings, which indicates the presence of QLR pair 
correlations. This extraordinarily difficult problem is tantamount to providing an analytical 
prediction of the MRJ state. Unfortunately, no such rigorous theory currently exists, even 



for the considerably simpler problem of predicting the scalar MRJ density |26[. The presence 
of QLR correlations makes this problem inherently nonlocal, and, therefore, methods that 
attempt to predict the MRJ state based only on packing fraction and local criteria, such as 



26). 



nearest-neighbor and Voronoi statistics, are invariably incomplete 

We have presented arguments in a recent letter to suggest strongly that hyperuniformity 
and quasi-long-range pair correlations are signatures of saturated MRJ packings of hard 
particles, including binary disks, ellipses, and superdisks j^. In this paper, we provide 
detailed evidence to show that polydisperse MRJ packings of hard spheres, though inho- 
mogeneous with respect to the number variance, possess local-volume-fraction fluctuations 
decaying faster than the volume of an observation window. This observation is consistent 
with the generalized Torquato-Stillinger conjecture that all strictly jammed saturated pack- 



ings of spheres are hyperuniform with respect to fluctuations in the local volume fraction. 
Our major results include: 



1. Infinite wavelength local number density fluctuations do not vanish for MR J pack- 
ings of polydisperse hard spheres. Local-volume-fraction fluctuations provide the ap- 
propriate structural description of these packings because they account correctly for 
the size distribution of the particles. Importantly, our work studying local-volume- 
fraction fluctuations contains previously-published results for 3D monodisperse MRJ 
hard sphere packings special case. 

2. Signature QLR pair correlations scaling asymptotically as r~( d+1 ^ in d Euclidean di- 
mensions are observed for all systems that we study, including binary disks with vary- 
ing size ratios and compositions and polydisperse disk packings. Our results suggest 
that these special correlations may be a universal feature of the MRJ state. 

3. Strict jamming places a strong constraint on the distribution of the available void 
space external to the particles such that hyperuniformity is observed when considering 
local-volume-fraction fluctuations, even when the point pattern of the sphere centers 
is locally inhomogeneous. 

4. The competition between maximal randomness and strict jamming of the packings 
ensures that the void-space distribution is sufficiently broad to induce QLR correlations 
between particles, thereby providing a direct qualitative structural explanation for the 
linear small-wavenumber region of the generalized scattering intensity. 



II. BACKGROUND AND DEFINITIONS 



A. Point patterns 



We consider point patterns to be realizations of stochastic point processes. Formally, a 
stochastic point process is a method of placing points in some space (such as M. d ) according to 
an underlying probability distribution. This random setting is quite general, incorporating 
cases in which the locations of the points are deterministically known, such as in a Bravais 



lattice 



28) 
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A statistically homogeneous point process is completely determined by the number den- 
sity p and the countably infinite set (in the thermodynamic limit) of n-particle correlation 
functions. The n-particle correlation function g n (ri, . . . , r n ) is proportional to the prob- 
ability density of finding n particles in volume elements around the positions r 1; ...,r n , 
regardless of the positions of the remaining particles in the system. For an arbitrary point 
process, deviations of g n from unity provide a measure of the correlations among points 
in the system. Note that specifying only a finite number M of the n-particle correlation 
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functions defines a class of microstructures with degenerate M-particle statistics 
particular interest is the pair correlation function g 2 , which defines the average number of 
particles surrounding a reference particle of the point process. Closely related to the pair 
correlation function is the total correlation function h(r) = ^( r ) — 1- Since g 2 (r) — > 1 
as r — > +oo (r = ||r||) for isotropic, translationally invariant systems without long-range 
order, it follows that ft(r) 4 in this limit, meaning that h is generally integrable with a 
well-defined Fourier transform. 

It is common in statistical mechanics when passing to reciprocal space to consider the 
associated structure factor S(k), which for a translationally invariant system is defined by 

S(k) = l + ph(k), (1) 

where h is the Fourier transform of the total correlation function, p is the number density, 
and k = ||k|| is the magnitude of the reciprocal variable to r. We utilize the following 
definition of the Fourier transform: 

/(k)= / /(r)ex P Hk-r)dr, (2) 

where kr = Yli=i ^i r i ^ s the conventional Euclidean inner product of two real- valued vectors. 
For radially-symmetric functions [i.e., /(r) = /(||r||) = f(r)}, the Fourier transform may be 
written 

B. Two-phase random heterogeneous media 

Closely related to the notion of a stochastic point process is that of a two-phase random 
heterogeneous medium (or random set), which we define to be a domain of space V C M. d of 



volume V < +00 that is composed of two regions: the phase 1 region Vi of volume fraction 



>i and the phase 2 region V2 of volume fraction 02 



30(. The statistical properties of each 



phase of the system are specified 
Sn\ which are defined by 



3V the countably infinite set of n-point probability functions 



|33j 




^)( ri ,...,r n ) = ( 
where 1^ is the indicator function for phase i 

/W(x)= r (5) 




The function S n defines the probability of finding n points at positions ri, . . . , r n all within 
the same phase. 

Upon subtracting the long-range behavior from S2, one obtains the autocovariance func- 
tion x(r) = ^(r) — 2 , which is generally integrable. It is important to recognize that the 
autocovariance function is independent of the choice of reference phase, meaning that it is 
a global descriptor of correlations within the system. This property will play a particularly 
important role in this paper when we consider the relationship between local-volume-fraction 
fluctuations and the void space between inclusions in the microstructure. The analog of the 
structure factor in this context is the so-called spectral density, which is the Fourier trans- 



orm x °f the autocovariance function 34j. The autocovariance function obeys the bounds 



[ul 



- min{(l - 0) 2 , 2 } < X (r) < (1 - 0)0, (6) 

where is the volume fraction of an arbitrary reference phase. We remark that it is an open 
problem to identify additional necessary and sufficient conditions that the autocovariance 



function must satisfy in order to correspond to a binary stochastic process 



Q0& 



C. Sphere packings and a categorization of jamming 

A sphere packing is a collection of non-overlapping spheres in <i-dimensional Euclidean 
space M. d . The packing density (equivalent to the volume fraction of the particle phase) 
is defined as the fraction of space covered by the spheres, which may be polydisperse. An 



S 



important characteristic of a packing is its degree of randomness (or the antithesis, order), 
which reflects nontrivial information of the packing structure. The degree of randomness 
(order) can be quantified by a set of order metrics ifj It is an open and challenging prob- 



lem to identify good order metrics, but it has recejat 
is itself a measure of order on large length scales 

One method of classifying sphere packings involves characterizing the extent to which 



been proposed that hyperuniformity 
61. 



particles are jammed. Torquato and Stillinger [39] , |40| have provided a precise definition of 
the term jamming and have proposed a hierarchical classification scheme for sphere pack- 
ings by invoking the notions of local, collective, and strict jamming. A packing is locally 
jammed if no particle in the system can be translated while fixing the positions of all other 
particles. A collectively jammed packing is locally jammed such that no subset of spheres 
can simultaneously be continuously displaced without moving its members out of contact 
both with one another and with the remainder set. A packing is strictly jammed if it is 
collectively jammed and if all globally uniform volume nonincreasing deformations of the 
system boundary are disallowed by the impenetrability constraints. The reader is referred 
to Ref. 



for further details. 



As previously mentioned, the maximally random jammed (MRJ) state is defined as the 
most disordered jammed packing in a given jamming category (i.e., locally, collectively, or 
strictly jammed) [lj. The MRJ state is well-defined for a given jamming category and choice 
of order metric, and it has recently supplanted the ill-defined random close packed (RCP) 
state [lj. In this paper, we focus on maximally random strictly jammed polydisperse sphere 
packings in M. d . 



D. Hyperuniformity in point processes: local number density fluctuations 

A hyperuniform point process has the property that the variance in the number of points 
in an observation window Q grows more slowly than the volume of that window. In the case 
of a spherical observation window, this definition implies that the local number variance 
a 2 N (R) grows more slowly than R d in d dimensions, where R is the radius of the observation 
window. Torquato and Stillinger have provided an exact expression for the local number 
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variance of a statistically homogeneous point process in a spherical observation window 



a 2 N (R) = pv{R) 



l + p h(r)a(r;R)dr 



(7) 



where R is the radius of the observation window, v(R) is the volume of the window, and 
a(r; R) is the so-called scaled intersection volume. The latter quantity is geometrically 
defined as the volume of space occupied by the intersection of two spheres of radius R 
separated by a distance r and normalized by the volume of a sphere v(R). Exact expressions 
for a(r; R) in arbitrary dimensions have been given by Torquato and Stillinger 41]. 

It is convenient to introduce a dimensionless density 4>, which need not correspond to the 
volume fraction, according to 

p7T d/2 D d 



= pv(D/2) 



(8) 



2 d r(l +d/2)' 

where D is a characteristic length scale of the system (e.g., the mean nearest-neighbor 
distance between points). The number variance admits the following asymptotic scaling |||: 



a%(R) 



A 



N 



D 



N 



d-l 



d-1 



(9) 



where o(x) denotes all terms of order less than x. This result is valid for all periodic point 
patterns (including lattices), quasicrystals that possess Bragg peaks, and disordered systems 
in which the pair correlation function $2 decays to unity exponentially fast 



for the asymptotic coefficients and are given by [5| 

l+p I h(r)dr = lim S(k) 

jRd ||k||->0 

pn(d) 



A N 
B N 



D 



h{v)\\v\\dT, 



5| . Explicit forms 

(10) 
(11) 



where n{d) = T(l + d/2)/{n 1 / 2 T[{d + l)/2]}. 

Any system with = satisfies the requirements for hyperuniformity. Although the 
expansion fl9]) will hold for all periodic and quasiperiodic point patterns with Bragg peaks, 
this behavior is not generally true for disordered hyperuniform systems. For example, it is 
known that if the total correlation function h ~ r -i d + 1 ) f or l ar g e r [S(k) ~ k for small k], 
then a N (R) ~ (a \nR + a^R^ 1 Such behavior occurs in maximally random jammed 
monodisperse sphere packings in three dimensions J4] and noninteracting spin-polarized 
fermion ground states Q]. Other examples of "anomalous" local density fluctuations 
have been characterized by Zachary and Torquato (]]. 
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E. Hyperuniformity in two-phase random heterogeneous media: local-volume- 
fraction fluctuations 



In order to define hyperuniformity for heterogeneous media we introduce the local volume 
fraction Tj(x) of phase i according to 

1 



T;(x; R) 



v(R) 



/ w (z)m;(z-x; R)dz, 



(12) 



where v (R) is the volume of the observation window and w is the corresponding indicator 
function. Using this definition, the variance o~%{R) in the local volume fraction is given by 



o*(R) 



v(R) 



x(r)a(r;R)dr, 



:i3i 



which is independent of the choice of reference phase. The variance in the local volume 



fraction admits the asymptotic expansion 

d 



P 



6] 



d+l 



d+l 



K = / X(r)dr = lim x(k) 

jRd l|k|K0 

K(d) 



D 



|r||x(r)dr. 



(14) 

(15) 
(16) 



The coefficients A T and B T in ( fl5i) and ( TlBT) control the asymptotic scaling of the fluctuations 
in the local volume fraction. It then follows that a 2 T decays faster than R~ d as R — >■ +oo for 
those systems such that 

lim_x(k) = 0, (17) 



IlklKO 



which defines a hyperuniform two-phase random heterogeneous medium |6] 



F. The effect of polydispersity on local fluctuations 



Here we consider how the presence of polydispersity in a sphere packing affects the fluctu- 
ations in the local number density and local volume fraction. For a packing of polydisperse 
spheres, the distribution of sphere radii is determined by a probability density f(R), and the 
fundamental statistical descriptors of the medium therefore involve averages 42| over 



the distribution of sphere sizes 



11] . For a packing of polydisperse spheres with M distinct 



11 



radii, the density function f(R) takes the form 



M 



f(R) = ^2 li 5(R-R i ), (18) 



i=l 



where 7, = Ni/N is the mole fraction of species i. As an example, we will continually refer 
in this paper to the average particle diameter (D)r of a polydisperse hard-disk packing: 



Al 



(D) R = J2l*D l , (19) 



i=i 



where Di is the diameter of species i. 

The autocovariance function y(r) for a heterogeneous medium consisting of impenetrable 
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polydisperse spheres is given by 

X(r) = p(f in t(r; + P 2 y (/i(x; #1, i? 2 )fint(r - x; JSi, R-i)) RlM dx, (20) 
which implies that 

A r = J X (r)dr (21) 
= p(v 2 (R)) R + p 2 (v{R 1 )v(R 2 ) [ h(r; R u R 2 )dr) . (22) 



Unlike for monodisperse sphere packings 



6|, the result (I22p shows that it is generally not 
possible to separate the shape information of the inclusions from the details of the point 
pattern generated by the sphere centers. This observation suggests that for polydisperse 
microstructures, hyperuniformity of the underlying point pattern does not induce hyperuni- 
formity with respect to local-volume-fraction fluctuations; conversely, it is also possible to 
find heterogeneous media for which A r = but An > 0. 

We note that the first term contributing to x{ r ) m (120]) can be interpreted as the proba- 
bility of finding two points, separated by a displacement r, in a single particle of the packing. 
The second term is therefore related to the probability of finding the points in two different 
particles. It is clear that only this latter term, containing the pair correlation function, can 
contribute to the linear small-wavenumber region of the spectral density, and it is therefore 
responsible, albeit in a highly nontrivial way, for the onset of QLR correlations in MRJ 
packings. 
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III. LOCAL DENSITY AND VOLUME FRACTION FLUCTUATIONS IN POLY- 
DISPERSE MRJ PACKINGS 



A. Generation of MRJ polydisperse hard disk packings 

Motivated by the generalized conjecture that all strictly jammed packings of d-dimensional 
spheres are hyperuniform with respect to local-volume-fraction fluctuations, we have gener- 
ated several packings of binary disks (2D) near the MRJ state using a modified Lubachevsky- 
Stillinger (LS) packing algorithm 45M49| . wherein particles with a fixed size ratio undergo 
event-driven molecular dynamics while simultaneously increasing in size according to some 
prescribed growth rate. The initial growth rate used in our simulations is 7 = 0.01, but near 
the jamming point, a much smaller growth rate 7 = 10~ 6 is used to establish well-defined 
interparticle contacts. The Lubachevsky-Stillinger algorithm has been shown to produce 
MRJ packings with these parameters consistent with the positively-correlated translational 



and orientational order metrics 



50 



5l| . Our statistics for binary packings are averaged 
over 50 configurations of 10000 particles; for polydisperse packings we have averaged over 
10 configurations of 10000 particles. Our results have been compared to systems of up to 
10 6 particles to verify invariance of the statistics with respect to system size. 

For this study, we have chosen a particle size ratio /3 = -Riargc/-Rsmaii = 1.4 and mole 
fractions 7 sm aii = 0.75 and 7i a r g o = 0.25; however, the focus of this study is to elucidate 
a universal property of polydisperse MRJ packings, and our results are expected to apply 
for a range of size ratios, size distributions, and mole fractions (see Figure [2]) 521 ] . For the 
remainder of this paper, we therefore focus on the binary case with the disclaimer that our 
results are expected to apply for general polydisperse MRJ packings. This algorithm has 



been used in previous studies of 3D monodisperse sphere packings near the MRJ state [47 ]. 
The resulting packings in this study have a final volume fraction ps 0.8475, which is below 
the the close-packed density cp = a/37t/6. Figure Q] (Section I) provides a typical realization 
of a binary packing; note that the particle distribution is saturated and nonperiodic. 

B. Hyperuniformity and local- volume- fraction fluctuations 

We have calculated both the structure factor S(k) and the spectral density x(k) for the 
binary MRJ disk packings using discrete Fourier transforms of the local density p(r) and 
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indicator function J(r) of the particle phase of the packing. Specifically, 



EjLi exp (~ik • ,0 
5(k) = — (k 0) 



Y$=i ex P ( _ik • r i) ™( k ; Rj) 



2 



X(k) = — (k ^ 0), (24) 



where 



m(k;R)= exp(-ik • r)Q(R - \\r\\)dr (25) 

is the Fourier transform of the indicator function for a d- dimensional sphere of radius R 
53I . Note that the shape of the enclosure (defined by a set of basis vectors {ej}) restricts 
the wavevectors such that k • e$ = 27m for all i, where n£Z. Since the zero wavevector is 
removed from the spectrum in the expressions above, one must define 

X(0)= lim x(||k||=A; min ), (26) 

N,V— H-oo 

with a similar expression for S(0). The limit here is taken at constant number density 
p = N/V, and /c m i n is the smallest computable wavevector as determined by the shape 
of the boundary. To obtain radially-symmetric forms of the structure factor and spectral 
density, we radially average over all wavevectors with equal magnitude. 

Our results for the structure factors and spectral densities of the binary and polydisperse 
MRJ hard disk packings are shown in Figure [2j The structure factors for these systems 
lacks Bragg peaks, reflecting the absence of long-range order. Of particular importance 
is the behavior of the structure factor near the origin. A fit of the small-fc region (k < 
0.5) with a third-order polynomial suggests S(0) ~ 0.104 > 0, meaning that the point 
pattern generated by the disk centers does not possess vanishing infinite-wavelength local 
number density fluctuations, an observation verified by direct calculation of the number 
variance in Figure [3j Note that the local number variance asymptotically scales with the 
volume of the observation window. This behavior differs from MRJ packings of monodisperse 

n 

spheres in three dimensions |47[ , where the structure factor decays linearly to zero for small 
wavenumbers. Bidisperse distributions of disks are inherently inhomogeneous with respect to 
the locations of the particle centers. For the systems studied here, the large concentration 
of small particles results in local clusters that are essentially close-packed. However, the 
introduction of large particles into the system generates effective grain boundaries between 
these local clusters, breaking the uniformity of the underlying point pattern. 
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FIG. 2. (a) Structure factor S(k) for a binary packing of disks near the MRJ state, (b) The 
corresponding spectral density x(k). Note that infinite- wavelength local density fluctuations are 
not suppressed unlike local volume fraction fluctuations on the equivalent length scale, (c) Spec- 
tral density for particle concentrations 7smaii = 

2/3 and 

Tiarge — 1/3 with size ratio f3 — 2.5. 
Hyper uniformity of the packing is unaffected by changes in these parameters as expected by the 
Torquato-Stillinger conjecture, (d) Spectral density for a polydisperse MRJ packing of disks with 
a uniform distribution of radii in the interval [R m i n , Rmax\- 



We note that the structure factor for our binary packing can be decomposed as 



S(k) = S s (k) + S L (k) + S c (k), 
15 



(27) 
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(a) (b) 
FIG. 3. (Color online) (a) Asymptotic fluctuations in the local number density for the MRJ binary 
disk packing, (b) Asymptotic fluctuations in the local volume fraction for the binary MRJ packing, 
demonstrating that the system is indeed hyperuniform. 



where 



E5i exp(-ik • x,- 



Ss(k) 
S L (k) 
S c (k) = 2Re 



N 



N 



E j =i exp(-ik • Tj) E£ ex P (* k ■ 



N 



(28) 
(29) 
(30) 



are the particle structure factors incorporating small- small, large-large, and small-large cor- 
relations, respectively. Note that N$ is the number of small particles with positions {r^}, 
and Nl is the number of large particles with positions {x^}. These partial contributions to 
the structure factor are shown in Figure HI None of the partial contributions to the structure 
factor possess a vanishing small-wavenumber region, implying that the local density fluctu- 
ations of the small and large particles each scale with the volume of an observation window 
as does the covariance. These observations clearly demonstrate that information contained 
in the structure factor is not sufficient to characterize the packings because it neglects the 
details of the particle shapes. 
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FIG. 4. (Color online) Partial contributions Ss(k), 5l(A;), and Sc(k) to the structure factor of the 
binary MRJ disk packings. 

A direct calculation of the spectral density for the binary MRJ system shows markedly 
different behavior from the structure factor as seen in Figure El One notices that infinite- 
wavelength local-volume-fraction fluctuations are suppressed by the system. We have fit the 
small- k region of the spectral density (in units of (-D)fj) using a third-order polynomial of the 
form ao + a\K + (I2K 2 + a^K 3 and have found ao = (1.0 ± 0.2) x 10~ 5 , strongly suggesting 
that this system is hyperuniform with respect to local-volume-fraction fluctuations. We 
have verified this claim by directly calculating the variance in the local volume fraction, 
shown in Figure |3j The local-volume-fraction fluctuations decay faster than the volume of 
an observation window and logarithmically slower than the surface area, consistent with the 
presence of QLR correlations. 

C. Properties of the spectral density of binary MRJ hard disk packings 

The spectral density is dominated by a peak with wavelength just below (D)r, approxi- 
mately corresponding to the small-particle diameter. This observation is consistent with the 
microstructure of the medium, which contains several regions of almost close-packed clusters 
of small particles. In order to further understand the behavior of the spectral density, it is 
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FIG. 5. The average magnitude {\ih(k; R)\)r of the spatial contribution to the spectral density. 

important to recognize that both the local distribution of inclusions and the shape informa- 
tion of the particles contribute to X- However, the relative influence of each component on 
the spectral density varies throughout the spectrum. To characterize this effect, we recall 
the representation (124jl of the spectral density, which admits the following upper bound: 



Figure [5] plots the average (\rh(k; R)\) R , which according to f l3~TT) controls the upper bound 
on the spectral density. Note that for small wavenumbers [k(D) r/ '(2tt) < 2] the bound 
( l3TT) places only weak constraints on the spectral density, implying that this region of the 
spectrum is controlled by information in the local structure of the heterogeneous medium. 
We emphasize that the shape information of the particles must still be included at small 
wavenumbers to account for vanishing infinite- wavelength local-volume- fraction fluctuations, 
but, as we show below, such geometric information only provides the appropriate weights for 
the partial contributions to the structure factor to induce hyperuniformity. In constrast, the 
large-wavenumber region of the spectrum closely follows the behavior of the upper bound 
( 13TT) . meaning that the length scale imposed by the decoration of the particle centers, here 
chosen to be disks, controls the spectral density. We emphasize that this portion of the spec- 
trum arises from the shape information of the inclusions themselves and is almost entirely 
independent of the distribution of particles in the system. Specifically, the particle indica- 



X(k) < Np {(\rh(k; R)\)r) 2 ■ 



(31) 
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FIG. 6. (Color online) Spectral density and corresponding shape contribution for a microstructure 
generated by the step-function g<i process. Note that the shape contribution controls the spectral 
density for large wavenumbers. 

tor function m(r;R), which has compact support [0,-R], possesses a Fourier representation 
rh(k; R) that is both long-ranged and has an intrinsic period, which for a binary packing 
leads to interference effects in the large-wavenumber region. 

Similar behavior arises in monodisperse systems, where the spectral density is exactly 
given x{k) = prh 2 (k; R) S (k) , and one can rigorously separate information contained in the 
point pattern generated by the sphere centers from the shape information of the inclusions. 
Figure |6] shows the spectral density and shape contribution for a system of impenetrable 
disks with pair correlation function g2{r) = Q(D — r). Note that the bound ( |3Tl) only 
applies when S(k) is at a maximum; for clarity we have omitted the associated scaling 
factor on the shape contribution. As with the binary MRJ packing, this "step-function 
process" possesses a spectral density that for small wavenumbers depends signficantly on 
the distribution of sphere centers; the large-wavenumber region is almost exactly equal to 
the shape contribution since the structure factor approaches an asymptotic value of unity. 
Although for binary systems polydispersity precludes the direct separation of the shape 
information from the underlying point pattern, these examples show that qualitatively these 
components continue to affect only specific portions of the spectral density. 
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x = R/<D> 

FIG. 7. (Color online) Standard deviation of the local volume fraction for the close-packed binary 
square and square lattices with the corresponding result for the MRJ binary packing. 

The dominance of the shape information for large wavenumbers has the surprising effect 
of almost completely suppressing fluctuations for K = k(D)n/ (2tt) ~ 2.5, thereby suggesting 
that local volume fraction fluctuations essentially vanish on this length scale. In actuality, the 
nonuniformity of the microstructure precludes a complete extinction of the variance in the 
local volume fraction; nevertheless, the local- volume-fraction variance (Figure [7]) undergoes 
a sharp change in slope near this length scale, highlighting a transition to the "geometrically- 
controlled" region of the spectral density. Physically, this behavior implies that observation 
windows with radii given by the appropriate wavelength capture the effective pore size 
surrounding the inclusions in such a way that the medium is essentially homogeneous on 
this local scale. This observation suggests that the void space surrounding the particles 
plays a central role in determining the spatial statistics of a heterogeneous system. 

D. Probing the origin of QLR pair correlations 

It is important to note that the small-A; behavior of the spectral density for the binary 
MRJ packing results in an "anomalous" asymptotic scaling of the variance in the local 
volume fraction. Specifically, our results suggest that the spectral density is nonanalytic at 
the origin with an expansion \(k) ~ aik + 0{k 2 ) as k — > 0, implying that the autocovariance 
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(a) (b) 

FIG. 8. (Color online) (a) Unit cell for the square lattice at close-packing, (b) Unit cell for a 
close-packed binary variant of the square lattice. 



function exhibits quasi-long-range behavior and scales with r~ 



in d dimensions]. 



Equivalent behavior has been observed for the structure factor of MRJ monodisperse sphere 
packings in three dimensions {4], and the number variance in that case has been shown to 
scale according to: 

a 2 N (R) ~ (&o + &ilni*)fl 2 + 0(fl). (32) 

A related scaling must also hold for the variance in the local volume fraction of our binary 
MRJ packings; specifically: 



<fla)~ fa+ ^' nJ?) + o(jr'). 



(33) 



One can directly see the effect of this behavior in Figure [7J where the variance in the lo- 
cal volume fraction for the MRJ packings is compared to the square lattice packing and 
its binary variant, the unit cells for which are shown in Figure [SJ Note that local volume 
fraction fluctuations in the MRJ packings decay more slowly than in either of the periodic 
systems; indeed, as the size distribution of the void space becomes smaller and more uniform, 
local-volume-fraction fluctuations are more rapidly suppressed on the global scale of the mi- 
crostructure. Since the integer packing is not even strictly jammed, this comparison suggests 
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that strict jamming is neither a necessary nor a strong determinant of hyperuniformity in 
heterogeneous media. 

Identifying the origin of the linear small-wavenumber region of the spectral density is 
an open problem that must be related to the structural features of the MRJ state. This 
problem is particularly difficult for our polydisperse packings since the underlying point 
pattern generated by the particle centroids possesses non- vanishing infinite-wavelength local 
number density fluctuations. However, the spectral density x(k) can be expressed in terms 
of the partial structure factors as 

X{k) = pm 2 (k; Rs)Ss(k) + pm 2 (k; R L )S L {k) + pm(k; R s )m(k; R L )S c {k), (34) 

where Rs and Rl are the radii of the small and large particles, respectively. This result 
follows directly from an expansion of fl2l]). As k — > 0, we therefore find 

X(0) = pv 2 {R s )S s {0) + pv 2 (R L )S L (0) + pv(R s )v(R L )S c (0), (35) 

implying that the particle volumes provide the appropriate weights to properly balance the 
small- and large-particle variances with the covariance between the particles. Additionally, 
since rh(k; R) possesses no linear term in its small-wavenumber Taylor expansion, it follows 
that the appearance a linear small-wavenumber region in the spectral density, and therefore 
QLR correlations, must involve an appropriate superposition of linear contributions in the 
partial structure factors. 

This observation suggests that quasi-long-range correlations in MRJ packings arise from 
the competing effects of strict jamming and maximal randomness. Incompressibility of the 
structure from strict jamming implies that particles should be correlated over several char- 
acteristic length scales; indeed, along the equilibrium branch of the binary hard-disk phase 
diagram, one expects that correlations become fully long-ranged at the close-packed density, 



corresponding to phase-separated lattice structures 54j. However, maximal randomness 
interferes with this long-range order, resulting in the apparent r~( d+1 ) asymptotic scaling 
in the pair correlation function. For our binary packings, this scaling is encoded in the 
partial pair correlations of the structure, which, after appropriate weighting with the shape 
information of the particles, induces hyperuniformity. 

We emphasize that a complete explanation for the appearance of the linear small- 
wavenumber region of the spectral density is intractable because the problem is inherently 
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FIG. 9. (Color online) The quadratic form (|35p . equal to the small- wavenumber limit of the spectral 
density x(0), as a function of the radii of the particles in the packing. A uniform decrease in the 
particle radii defines a basin of hyperuniform systems. 



non-local due to hyperuniformity and the presence of QLR correlations. Furthermore, it has 
recently been established that decreasing the exponential scaling of the small-wavenumber 
is associated with greater disorder within a many-particle system, potentially even induc- 
ing clustering among the particles {55}], an effect clearly inconsistent with strict jamming. 
More generally, hyperuniformity is associated with an effective interparticle repulsion that 
attempts to evenly distribute the particles throughout space, and the length scale of this 
repulsion increases with increasing exponential scaling of the small-wavenumber region of 



the spectral density [55|. The conditions of saturation and impenetrability in an MRJ 
packing induce structural order that apparently competes with the constraint of maximal 
randomness to minimize the scaling of the small-wavenumber region of the spectral density 
to its smallest integer value, an effect that we argue is physically tied to the void-space 
distribution of the system. 

To motivate our discussion of the void space and its fundamental importance to polydis- 
perse MRJ packings, we mention a peculiar property of hyperuniformity in these packings. 
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Namely, one can maintain hyper uniformity in polydisperse MRJ hard-particle systems even 
upon shrinking the particles at a fixed size ratio so long as one does not affect the under- 
lying statistical distribution of the point process; this effect is illustrated in Figure [EH To 
understand this behavior, we note that the leading-order term governing the expansion of 
the Fourier transform of the particle indicator function rh(k; R) is the volume of the particle 
ttR 2 , meaning that for small wavenumbers we have (in the case of the binary packings) 

2 

(36) 



X(k) ~ - 



Nsmall Marge 



j=l l=\ 

which, upon rescaling Ri — > kRi for re < 1, suggests 

Xshrink(k) ~ re 4 XMRj(k) (k 0). (37) 

Therefore, hyperuniformity of the MRJ binary packing is not lost when performing this scal- 
ing operation. Physically, the small-wavenumber region of the spectral density effectively 
homogenizes the medium due to the coupling between the wavenumber and the particle 
radius in the indicator function rh(k; R), meaning that this region is not affected by changes 
at the boundaries of the particles so long as the underlying statistics of the point pattern 
generating the medium remain constant. This observation suggests that a saturated and 
strictly jammed sphere packing extends naturally to an uncountably infinite family of hy- 
peruniform heterogeneous media related by an appropriate scaling parameter. One can 
map this behavior directly using the result (135]) . which indicates that the spectral density 
at small wavenumbers defines a quadratic form in the particle radii when the underlying 
point pattern of the disk centers is held fixed; the corresponding curve is plotted in Figure 
[9j Furthermore, this scaling has the effect of deforming the void space surrounding the 
particles in a uniform manner that preserves the regularity of the pore distribution even 
upon breaking both the jamming and saturation constraints. The medium therefore is able 
to retain information from the strictly jammed configuration even upon relaxing the sizes 
of the inclusions. Additionally, it follows from the quadratic form ( 1351) that non-uniform 
changes in the particle radii inherently break the hyperuniformity of the packing. 

This behavior has important implications for MRJ packings. First, this effect is geometric 
in origin and depends on the explicit inclusion of particle shape information to appropri- 
ately balance the partial structure factors of the packing. This observation immediately 
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FIG. 10. (Color online) Spectral densities upon shrinking the particles in the binary MRJ packings 
at fixed size ratio /3 = i?iaxge/^small- The parameter k = i? s hrink/-^MRj is the scaling factor. 



implies that "point" information contained in the particle centers is not sufficient to de- 
scribe the system. Second, the presence of interparticle contacts is not essential for the 
onset of hyperuniformity and QLR correlations. This subtle point, though easy to appre- 
ciate mathematically, is highly nontrivial from a physical perspective. It implies that on 
the global scale of the microstructure, the space is regularized in such a way that hype- 
runiformity is preserved upon making uniform scaling deformations of the particle phase. 
Importantly, since the spectral density is a descriptor of both the particle and the void 
phases, this homogenization is invariant to the choice of reference phase. Note in particular 
that upon breaking the contact network of the MRJ packing, the void space can become 
connected throughout the space, implying that the spatial statistics of the microstructure 
must account for correlations over large length scales. It is for these reasons that we turn 
our attention to the void-space distribution in the following sections. 

Note that our scaling analysis does not apply if the particle radii are uniformly increased 
since breaking the impenetrability constraint implies that higher-order microstructural in- 
formation is necessary to characterize the medium, and ( 1361) is no longer valid. Additionally, 
if the particles are shrunk to the point that the small particle radius vanishes, then the 
packing will no longer be hyperuniform since important shape information about the sys- 
tem has been lost. In terms of the void space, the distribution of pore sizes will be skewed 
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toward higher values due to the sudden appearance of "holes" in the microstructure, corre- 
sponding to the lost small particles, thereby de-regularizing the microstructure and breaking 
hyp erunif ormity. 



IV. CHARACTERIZATION OF THE PORE-SIZE DISTRIBUTION 

A. Numerical evaluation of pore-size statistics 

It is clear from the discussion above that the void phase surrounding the disk inclusions 
of the binary MRJ packings plays a significant role in characterizing local fluctuations in the 
medium. Indeed, we argue that the conditions of strict jamming and saturation place strong 
constraints on the distribution of the pore sizes, effectively regularizing the local structure 
around each disk such that the system is hyperuniform with respect to local-volume-fraction 
fluctuations even though the disk centers do not consitute a hyperuniform point pattern. To 
support our arguments concerning the void space, we have quantified the size of the available 
void space using the so-called complementary pore-size cumulative distribution function F(5) 



111 ], which represents the fraction of the void space external to the inclusion phase with a 
pore radius larger than 5. Equivalently, if we define P (8) as the probability density that a 
randomly chosen point in the void space lies within 5 and 5 + dS from the nearest point on 
the void-inclusion interface, then: 

r+oo 

F{5) = / P(A)dA. (38) 



J s 

Figure (TT] shows the cumulative pore-size distribution function for the binary MRJ pack- 
ings along with the corresponding results for the integer and binary integer packings and a 
system of equilibrium monodisperse hard disks at volume fraction <fi = 0.5. The equilibrium 
hard disk packing is neither saturated nor jammed; it is also known that this system is 
not hyperuniform jlj]. As a result, it possesses a broad distribution of pore sizes, and the 
probability of finding large pores (i.e., on the order of a particle diameter) is nonvanish- 
ing. In contrast, both the integer packing and its binary variant possess narrow pore-size 
distributions with compact supports, which result from the close-packed nature of the pack- 
ings. Interestingly, Figure [11] indicates that the binary MRJ packing also possesses a narrow 
pore-size distribution that essentially falls between the square and binary square packings. 
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FIG. 11. (Color online) Cumulative pore-size distributions F for the binary MRJ packing, the 
square lattice Z 2 , a saturated binary variant of Z 2 , and a system of equilibrium hard spheres with 
volume fraction A = 0.5. 



This observation suggests that the void space is highly constrained by the condition of strict 
jamming and almost regular in its distribution. 

One can directly compute the moments (5 n ) of the pore size density from the cumulative 
distribution F(S) using the relation 111 ]: 



r+oo 

(S n ) =n 6 n - 1 F(S)dS. 
Jo 



(39) 



Lower-order moments of the pore-size density arise in bounds on the mean survival times 
and principal relaxation times of heterogeneous materials and are thus important descriptors 
of microstructures {|36| [stJ (see Section V.B below). Table [J provides the first moments and 
standard deviations of the pore size for the systems mentioned in Figure [TTJ The average 
pore sizes for the hyperuniform systems are significantly smaller than for the equilibrium 
hard disk packing. In fact, the pore size distribution of the binary MRJ packing appears 
to be more localized about its mean than even the integer packing. We conclude from this 
information that the available void space in the binary MRJ disk packing is sufficiently 
restricted such that local volume fraction fluctuations decay faster than the volume of an 
observation window, thereby reflecting the underlying regularity of the pore space. 
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TABLE I. First moments and standard deviations of the pore size for distributions of hard disks. 



(8) /(D) <rs/(D) 

binary MRJ 0.0438 0.0380 

1? 0.0626 0.0480 

binary 1? 0.0332 0.0258 

equil. HS (</> = 0.5) 0.1612 0.1267 



B. Bounds on the pore size distribution 



It is important to note that the pore-size distribution for a hyperuniform hard-particle 
packing with volume fraction must asymptotically decay faster than the corresponding 
distribution for a non-hyperuniform medium with the same o ne-p oint statistics; this claim 
is a generalization of an argument by Gabrielli and Torquato 15[. Since non-hyperuniform 
media are by definition "Poisson-like," it suffices then to consider the pore-size distribution 
for a polydisperse system of fully penetrable disks. The cumulative pore-size distribution 



Fp(S) for such a system is known analytically 



11 



l33| and is given by 



F P (S) = (l/0i)exp [-p{ v (S + R)) 



if 



(40) 



where the angular brackets denote an average over the disk radii and <f>\ is the volume 
fraction of the matrix phase. We emphasize in these calculations that we are considering 
the pore-size distribution for a 2D Pois son- distributed heterogeneous medium with volume 
fraction <j) that is the same as a reference hyperuniform system. However, equivalence of 
the one-point statistics does not imply that the number densities are the same between the 
systems. To account for this discrepancy, we write p = ri/{ir(R 2 ) R ) = ln(l / 4>i) / (R 2 ) r), 
which implies 



F P (6) = exp —(s{D) R /(j)i)5 2 - (s(D) R /^)6 



(41) 



where s = p<j)i(s(R)) R is the specific surface of a Poisson-distributed medium and s(R) is 
the surface area of a d- dimensional sphere of radius R. Note that we have introduced the 
length scale (D) R and reduced variable 5 = S/(D) R in (HTj) . Letting b = s(D) R /(j) 1 , we find 
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the following upper bounds for the first and second moments of the pore size density: 

(*>ul = exp [6/4] erfc (Vb/2) (42) 

(5 2 )ul = ^(i-^)ul)- (43) 

The above Poisson bound will be rigorously true for the full pore-size distribution of 
our polydisperse MRJ packings based solely on the presence of hyperuniformity. However, 
we can significantly tighten the bounds of the moments of the pore-size distribution for 
binary MRJ packings based on our analysis below of the local voids. Specifically, since the 
distribution of voids within the MRJ packings is dominated by three- and four-particle loops, 
the average pore size must be less than the corresponding average pore size of a saturated 
integer lattice. This effect is explicitly shown in Figure [HI where it is clear that the integer 
lattice provides a good approximate upper bound for the full pore-size distribution with the 
exception of the tail, due to the presence of higher-order loops in the MRJ packing. Since 
this tail must decrease faster than exponentially with a cut-off at the small-particle radius 
(because of saturation), we therefore have the following tighter upper bounds on the first 
and second moments of the pore size distribution: 

(5) < (5) & (44) 
(S 2 ) < (S 2 ) Z 2. (45) 

It is also possible to find a simple lower bound on the pore-size distribution that is 
applicable for arbitrary heterogeneous media (hyperuniform or not). Specifically, we utilize 
the following series representation for the cumulative pore-size distribution fll| : 



01 



(46) 

It is well-known that keeping terms up to order in (j4"6|) places an upper bound on F(5) 
for k even and enforces a lower bound for k odd Therefore, by expanding the series to 
order p we obtain the following lower bound on F(S): 

F(5) > max + fl)> fl ],oJ; (47) 

the max operation in (T47|) enforces the trivial lower bound F(S) > 0. Note that ( )47|) 
is equivalent to the O(p) expansion of the pore-size distribution for a Poisson-distributed 
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medium; however, the number density (and therefore the specific surface) of the system is 
not necessarily the same as the Poisson medium. We will now assume that our hyperuniform 
reference medium consists of impenetrable spheres with volume fraction = p(v(R)) R ; we 
also remark that the definition of the specific surface for impenetrable spheres s = p(s(R)) R 
is different than for a fully-penetrable system. It now follows that 

F(6) > max jl — b5 2 - bS, j , (48) 

where b = s{D) R /<f>\ as in the Poisson case. The lower bound (148]) first reaches zero at 



2 

We thus obtain the following lower bounds on the first and second moments of the pore size 
density 

(S>11 = (50) 

Figure [T2] plots the upper and lower bounds ( 142]) and (150]) on the mean pore size for hy- 
peruniform binary heterogeneous media at volume fraction 0. We emphasize that these 
bounds account only for the hyperuniformity of the packings and thus place constraints on 
the pore-size distributions of any binary hyperuniform hard-particle packing. 
For our binary MRJ packings, the specific surface is 

(7small + 7large/3) 2 



SMRJ<£V01 = (HMD) 2 R /(R 2 ) R = (40/0 ly 

. 7small + 7largeP 

The specific surface for the equivalent Poisson-distributed medium is 



21.6917. (52) 



s P {D) R /<^ = -H<f>i)(D) 2 R /(R 2 ) R « 7.34037; (53) 

as expected, this value is less than the corresponding result for the impenetrable case. Using 
these parameters, we collect in Table HT1 the bounds on the pore size of our binary MRJ disk 
packings. We immediately notice that the numerical values are well within the bounds 
above; it follows that the one-point statistics of the medium do not place strong constraints 
on the distribution of the void space. Furthermore, this evidence suggests hyperuniformity 
strongly regularizes the pore space compared to the Poisson-distributed system. 
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— Poisson upper bound 

— Series lower bound 




FIG. 12. (Color online) Upper and lower bounds on the mean pore size applicable for binary 
hyperuniform heterogeneous media at volume fraction <f> with the same composition as the MRJ 
binary disk packings studied here. 



TABLE II. Bounds on the moments of the pore-size distribution for binary MRJ packings of hard 
disks. 



(S) (5* 



Poisson upper bound 0.1135 0.0227 
I? upper bound 0.0626 6.219 xl0~ 3 
Series lower bound 0.0224 6.635xl0" 4 



MRJ binary packing 0.0438 3.362x10 



-3 



It is interesting to note that the series lower bound provides a better estimate of the 
mean pore size than the Poisson upper bound. This observation is reasonable since the 
compact support of the lower bound more closely matches the behavior of the actual pore- 
size distribution, which also has compact support due to the saturation of the packing. 
However, it is unclear if the constraint of saturation can be relaxed while still maintaining 
hyper uniformity in the medium. One can conjecture that so long as the probability of 
finding a large pore decays faster than the corresponding behavior for a Poisson pattern, 
then hyperuniformity holds. 
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Although our results provide a quantitative basis for understanding the appearance of 
hyper uniformity in polydisperse MRJ packings, a complete explanation for the linear scaling 
of the small- wavenumber region of the spectral density is still not apparent. Again, we 
emphasize that this problem is inherently non-local, and its solution must account for the 
presence of QLR correlations within the packing. Such correlations are difficult to discern 
with the distribution F(8) that we have presented here, which is essentially a one-point 



descriptor of heterogeneous media 



111- 



In the following section, we characterize the allowable n-particle loops within a binary 
MRJ packing. This analysis suggests that the variance in pore sizes is bounded by the strict 
jamming of the packings, supporting our argument these constraints are sufficient to induce 
hyp er uniformity. 



V. CHARACTERIZATION OF n-PARTICLE LOOPS IN MRJ BINARY DISK 
PACKINGS 

The preceding discussion shows how the void space is a fundamental descriptor of a two- 
phase random heterogeneous medium. As we have seen, fluctuations in the local volume 
fraction contain information about the void space and therefore provide a more complete 
picture of hyperuniformity in heterogeneous media. Here we quantitatively characterize the 
voids in MRJ binary disk packings and argue that strict jamming restricts the variance of 
pore sizes. 

An n-particle void is associated with a loop of n contacting disks (an n-particle loop) in 
which each disk only contacts two neighbors. These loops are defined topologically, meaning 
that their identification is invariant to local shears of the particle contacts. The smallest 
loops contain three mutually contacting particles. The constraint of saturation places an 
upper bound on the number of particles in a loop, and the largest loop we find in the MRJ 
binary disk packings contains 6 particles. The number of n-particle loops in a packing 
decreases rapidly as n increases. 

The area of the void associated with a loop can be rigorously computed by substracting 
the area of particles falling into the polygon constructed by connecting the centers of the 
disks in the loop. We note that, except for the 3-particle loop, the polygons associated 
with other n-loops (n > 4) can be continuously deformed while maintaining the contacts 
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FIG. 13. (Color online) The four distinct 3-particle loops in a binary circular disk packing. 

between the particles. Such deformations change the area of the voids associated with the 
loop. The rigidity of MRJ packings also strongly constrains the number of particles in a 
loop. In the MRJ packings of binary disks, the majority of the voids involve 3-particle 
loops, which form the "backbone" of the network. Certain 4-particle loops can be observed 
at the effective grain boundaries between disks with different sizes. A large portion of these 
4-particle loops are very distorted with a void area almost equal to that associated with 
two 3-particle loops. Loops with more particles are rare in the packings, and the observed 
ones are all strongly distorted with void areas almost equal to those associated with 3- or 
4-particle loops. Therefore, we only discuss the voids associated with 3- and 4-particle loops 
here, focusing on the maximum possible void areas associated with these loops. 

There are four distinct 3-particle loops in the binary disk packings as shown in Figure [T3l 
Defining the radii of small and large disks be R\ and i?2, respectively, the n-particle void 
areas X na , normalized by the volume of a small particle, are \ 3 a ~ 0.0513, X^b ~ 0.0631, 
A3C ~ 0.0790 and A3D ~ 0.1060 for a size ratio (3 = R2/R1 = 1.4. It can be seen clearly that 
the void areas constitute a small fraction of the small particle area. Since these 3-particle 
voids dominate the packing, it is reasonable to conclude that the regularity of these local 
clusters enforces hyperuniformity on the medium despite the nonuniform distribution of the 
sphere centers. 

There are total of six distinct types of 4-particle loops (see Figure UM in the binary disk 
packings, each associated with a void that can possess a spectrum of shapes and sizes by 
distorting the quadrangle formed by the centers of the disks. Although most of the void areas 
associated with these highly distorted 4-particle loops are almost equal to the area associated 
with two 3-particle loops, there are still relatively large voids. The maximum normalized void 
areas are given by \ 4A ~ 0.2732, \ 4B « 0.3208, A 4C « 0.3934, A 4 d ~ 0.3782, \ AE « 0.4498 
and A 4 i? ~ 0.5355. These values are relatively large but more tightly distributed around the 
mean when compared to the normalized areas for 3-particle voids. As a result, they can be 
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FIG. 14. (Color online) The six distinct families of 4-particle loops in a binary circular disk packing. 

considered as an estimate of the upper bound on the degree of local inhomogeneity that can 
be consistent with hyperuniformity of the packings. 

VI. CONCLUDING REMARKS 

We have provided a detailed study of local- volume- fraction fluctuations in MRJ packings 
of polydisperse hard disks and have shown that these systems are hyperuniform with quasi- 
long-range correlations. Our results strongly suggest that these QLR correlations are a 
signature of MRJ hard-particle packings, in contrast to previous misconceptions in the 



literature [24j. Although it is true that the structure factors for these systems do not vanish 
at small wavenumbers except in the monodisperse limit, we have shown that the more 
appropriate structural descriptor of MRJ packings is the spectral density, which accounts 
appropriately for the shape information of the particles. Our work therefore generalizes 
the Torquato-Stillinger conjecture for hyperuniform point patterns by suggesting that all 
saturated, strictly-jammed sphere packings are hyperuniform with respect to local- volume- 
fraction fluctuations. Furthermore, MRJ sphere packings are expected to exhibit quasi- 
long-range pair correlations scaling as r _ ^ +1 ^ in d Euclidean dimensions. Importantly, this 
generalization contains previously-published results for monodisperse MRJ packings 4j as 
the special case where the particle shape information can be rigorously separated from the 
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"point" information of the particle centers. 

Based on the observations that local-volume-fraction fluctuations are invariant to the 
choice of reference phase and that one can maintain hyperuniformity in MRJ polydisperse 
packings under a uniform scaling deformation of the particles, we have argued that the 
onset of hyperuniformity and quasi-long-range correlations results from the homogenization 
of the void space external to the particles. The conditions of saturation and strict jamming 
limit the sizes and shapes of the local voids, which are completely determined by the contact 
network of particles. Using a local analysis of the voids and rigorous bounds on the pore-size 
distribution, we have argued that the void space is highly constrained by strict jamming, 
thereby suppressing infinite-wavelength local-volume-fraction fluctuations. Furthermore, we 
suggest that the presence of quasi-long-range correlations reflects the inherent structural 
correlations from the contact network between void shapes. Specifically, saturation and strict 
jamming of the packings compete with the maximal randomness of the particle distributions 
to drive the small-wavenumber region of the spectral density to its smallest integer value. 

Although our work has addressed important problems related to the structural proper- 
ties of MRJ hard-particle packings, a number of unanswered questions remain. First, a 
rigorous foundation for the observed linear scaling in the small-wavenumber region of the 
spectral density is still missing. This problem is immensely difficult to handle theoretically 
because the problem is inherently long-ranged, and any local analysis of the MRJ structure 



will therefore be unable to account correctly for this behavior |26|. Second, although we 
have considered only polydisperse MRJ packings of <i-dimensional spheres, our results are 
expected to hold more generally for strictly jammed packings of hard particles of arbitrary 
geometry. Indeed, our arguments concerning the void space of an MRJ packing are easily ex- 
tended to include these more general cases, and in a companion paper we will provide direct 



evidence that MRJ packings of hard ellipses and superdisks [58[ are also hyperuniform with 
signature quasi-long-range correlations indicated by linear scaling of the small-wavenumber 
region of the spectral density. These results will suggest a remarkably strong extension 
of the Torquato-Stillinger conjecture, namely that all maximally random strictly jammed 
saturated packings of hard particles, including those with size- and shape-distributions, 
are hyperuniform with universal quasi-long-range correlations [via the two-point probability 
function S2 (r)] scaling asymptotically as r~^ d+1 \ 

We mention that recent, independent work, appearing in preprint form shortly after our 
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own manuscript was posted, [59J has shown that hyperuniformity and vanishing infinite- 
wavelength local-volume-fraction fluctuations in saturated MRJ packings of polydisperse 
disks are consistent with certain fluctuation-response relations involving a generalized "com- 
pressibility." However, this work does not address the appearance of quasi-long-range pair 
correlations in MRJ hard-particle packings and does not consider MRJ packings of non- 



spherical particles as in our companion paper 
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Our arguments suggest in particular that certain quantum many-body systems and cos- 
mological structures with the same linear small-wavenumber scaling in the structure factor 
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10] are statistically "rigid" in the sense that their microstructures are effectively homo- 
geneous over large length scales with vanishing infinite-wavelength number density fluctua- 
tions. These unique features are inherently linked to the structural properties of the system, 
independent of the physical model itself. 

Finally, we note that the systems examined in this paper are constrained to be both sat- 
urated and strictly jammed. Saturation of the packings is responsible for enforcing compact 
support in the pore-size distribution function and therefore plays an important role in regu- 
larizing the void space surrounding the jammed disks. However, it is not clear if saturation 
is a necessary condition to ensure hyperuniformity in strictly jammed heterogeneous media; 
this situation corresponds to a weakening of the Torquato-Stillinger conjecture. Namely, 
what conditions must a strictly jammed but unsaturated packing of hard spheres meet in 
order to be hyperuniform? We note that the event-driven molecular dynamics algorithm 
used here to generated the binary MRJ disk packings inherently precludes the presence of 
arbitrarily large "holes" within the packing. However, it does produce a small concentra- 
tion (~2.5%) of "rattler" particles, which are particles that are free to move within some 
small caged region of the packing. Removal of such particles is known to break hyperunifor- 
mity of the medium even though the strict jamming of the surrounding structure remains 
4j; therefore, strict jamming alone is not sufficient to induce hyperuniformity if large holes 
are "common enough" in the statistics of the microstructure. This scenario corresponds to 
skewing the pore-size distribution and thereby de-regularizing the void space. 

It is indeed possible to construct strictly jammed packings of two-dimensional disks with 
a hole of arbitrarily large size by starting with a ring of particles encompassing a large pore; 
this ring can be jammed by surrounding it with a close-packed collection of particles that 
approximates an impenetrable "wall" (see Figure [T5lh One is then free to construct any 
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FIG. 15. (Color online) A jammed configuration of particles (up to rotations) surrounding a circular 
void. By replacing the wall with a close-packed collection of particles, one can construct strictly 
jammed packings in the plane with pores of arbitrarily large size. 



strictly jammed system of particles outside the hole [60(. If such holes are sufficiently rare, 
then it is possible that the system may still be hyperuniform since the pores do not become 
very large on average. 
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